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5 FIELD OF THE INVENTION 

This invention relates to the design and analysis of interconnections between 
electrical systems. More specifically, the invention relates to the efficient modeling and 
simulation of multiconductor transmission lines by automatically selecting macromodels 
1 0 for any given transmission line structure. 

BACKGROUND OF THE INVENTION 

High-performance electrical system design requires the use of sophisticated 
interconnections between system components. These interconnections must be designed 
so as to achieve three interrelated objectives: minimize signaling delay between 
components; minimize electromagnetic cross talk noise between interconnections; and 
ensure immunity to external electromagnetic interferences. Two main examples of these 
interconnections are multiconductor transmission lines and connectors. The present 
invention is concerned mainly with multiconductor transmission lines. 

Multiconductor transmission lines are present throughout any electrical system 
comprising several integrated circuits (chips). They are used on a printed circuit board 
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(PCB) for signal transmission between different chipsets. Multiconductor transmission 
lines are also used to transfer signals inside a package containing the chip as well as to 
transfer signals inside the chip itself. PCB and package transmission lines are known in 
the art as off-chip transmission lines while the lines responsible for transferring signals 
within the chip are known in the art as on-chip transmission lines. The transmission line 
nature of an interconnection depends on the wavelength of the signal carried by the 
interconnection. With the constant increase in electronic system speed, the signal 
wavelengths are becoming shorter. The net result is that more and more of the 
interconnections are behaving as transmission lines, which makes the task of modeling 
and analyzing overwhelming. The situation is rendered even more complex by the fact 
that the on-chip and off-chip transmission lines behave very differently in terms of the 
losses (i.e., attenuation and the distortion) incurred by the signals they carry. Because of 
their small cross sections, on-chip transmission lines are very lossy relative to off-chip 
transmission lines. Among off-chip transmission lines, packaging transmission lines are 
in general more lossy than PCB transmission lines. These differences in location (on- 
chip, off-chip), length (short, long), losses (high, low), and signal wavelength make the 
efficient modeling, simulation, and analysis of transmission lines a difficult engineering 
task. 

There are two main macromodeling approaches used to analyze the behavior of a 
transmission line. The first approach is based on a preliminary extraction of the pure 
delay (also called time-of-flight delay) incurred by the signal as it is transmitted. An 
instance of such approach is the method of characteristics (MoC) (Branin, IEEE Proc. 
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Vol. 55, pp. 2012-2013, 1967) and its various generalizations (Gruodis and Chang, IBM 
J. of Res. Dev. Vol. 25, pp. 25-41, 1982). The second macromodeling approach 
represents the transmission line with a cascade of electrical cells, each cell comprising 
lumped circuit elements such as resistors, capacitors, and inductors. This approach is 
equivalent to approximating the transfer function of the transmission line, which is a 
transcendental function, with a rational (i.e., non transcendental) function (e.g. Dounavis 
et. al., IEEE Transaction on Advanced Packaging, Vol. 22, pp. 382-392, 2000). The first 
macromodeling approach is referred to as the delay extraction approach while the second 
macromodeling approach is referred to as the rational approximation approach. 

It is well known in the art (e.g., Elfadel et. al., IEEE Transactions on Advanced 
Packaging, Vol. 25, pp. 143-153, 2002) that long, low-loss lines, such as coaxial cable 
(e.g., those connecting processing nodes in a supercomputer) are efficiently simulated 
using macromodels based on delay extraction approach. Short, high-loss transmission 
lines, such as on-chip bus lines (e.g., global bus connecting cable and CPU in a 
microprocessor) are efficiently simulated using macromodels based on rational 
approximation approach. Commercial circuit simulators, such as HSPICE, offer both 
types of modeling approaches to users. However, those simulators require the user to 
select which approach to use in a given situation. Users lacking expertise in transmission 
line theory may select the wrong or less efficient transmission line model (e.g., rational 
approximation model for a long, lossless line), thus incurring a significant cost in terms 
of transmission line simulation efficiency and accuracy. There is therefore a need to 
develop an automatic selection system and method for transmission line macromodels 
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based solely on the physical and geometric characteristics of the transmission line. Such 
an automatic selection system can be part of a computer-aided-design (CAD) tool, such 
as a circuit simulator, that will handle all matters related to the efficient and accurate 
simulation of transmission lines throughout the electrical system. Such automatic 
5 selection will also be crucial for system-level simulation for high-performance electronic 
systems (e.g., mainframe computers) where the number of transmission lines is very large 
and it is not known a priori what type of macromodel is most appropriate for a given 
transmission line. The present invention satisfies these needs by providing a system and 
method for the automatic selection of transmission line macromodels. 

10 

SUMMARY OF THE INVENTION 

The automatic selection is based on the transmission line length, the line per-unit- 
length parameters, and the maximum frequency of operation. In accordance with the 
1 5 teachings of the present invention, a method and system are provided which analyzes the 
exponential matrix solution of the Telegrapher's equation of a transmission line in order 
to select the better and more efficient macromodel to be used for simulating 
multiconductor transmission lines. 

20 The analysis considers the lossy transmission line as a perturbation of the lossless 

transmission line, which behaves as an ideal delay for which a delay-based macromodel 
gives the exact answer very efficiently. Taylor series analysis is then used to derive a 
mathematical criterion from which one can decompose the behavior of the lossy line into 
a pure delay portion and a pure loss portion. This criteria is then used to derive a critical 
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line length defining the boundary between the two macromodel efficiency regions : 
rational approximation macromodel for line lengths below the critical length and delay 
extraction macromodel for lines of length above the critical length. 

An aspect of this invention is, therefore, the provision of an improved system and 
method for a computer electrical system design. Another aspect of the present invention 
is the provision of an improved system and method for computer interconnection system 
design. A further aspect of the present invention is the provision of an improved system 
and method for computer interconnection system design using macromodeling of 
multiconductor transmission lines. A still further aspect of the present invention is the 
provision of an improved system and method for computer interconnection system design 
using automatic macromodel selection of a multiconductor transmission line. A yet 
further aspect of the present invention is the provision of an improved system and method 
for automatic macromodel selection of a multiconductor transmission line that can 
incorporated in a computer-aided-design tool. 

This invention provides a system, a method and a signal bearing medium that 
tangibly embodies a program of machine-readable instructions executable by a digital 
processing apparatus to perform operations to automatically select a macromodel from a 
set of macromodels for use in simulating a transmission line. The operations include 
providing input parameters of per-unit-length resistance (R), per-unit-length inductance 
(L), per-unit-length conductance (G), per-unit-length capacitance (C), length (d), and 
maximum operating frequency (co ma ^; computing total distortion (A T ) from the input 
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parameters; providing an error threshold (e); comparing the total distortion with the error 
threshold; and automatically selecting a macromodel based upon whether the total 
distortion is more or less than the error threshold. ■ 

This invention further provides a system, a method and a signal bearing medium that 
tangibly embodies a program of machine-readable instructions executable by a digital 
processing apparatus to perform operations to automatically select a macromodel from a 
set of macromodels for use in simulating a transmission line. The operations include 
providing input parameters of per-unit-length resistance (R), per-unit-length capacitance 
(C), length (d), and maximum operating frequency (ow), and error threshold (e); 
computing a critical length (d crit i Ca i) from the input parameters; comparing the length of 
the transmission line with the critical length; and automatically selecting a macromodel 
based upon whether the length of the transmission line is less than or greater than the 
critical length. 

BRIEF DESCRIPTION OF THE FIGURES 

Figure 1 is a schematic diagram of a transmission line represented as a two-port 

circuit. 

Figure 2 is a schematic electrical circuit diagram representation of a typical cell of 
a. rational approximation macromodel of a transmission line. 
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Figure 3 is a schematic electrical circuit diagram representation of a typical cell of 
a delay extraction macromodel of a transmission line. 

Figure 4 is a graphical representation summarizing the CPU efficiency of the 
5 rational approximation macromodel and delay extraction macromodel of a transmission 
line. 

Figure 5 is a flow chart of the overall CAD process for simulation of an electrical 
system with several multiconductor transmission lines. 

10 

Figure 6 is a block diagram of a preferred embodiment of an automatic 
macromodel selection process incorporating the present invention used in the CAD 
system in Figure 5. 

15 Figure 7 is a flow chart of an automatic transmission line macromodel selection 

process comprising the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

20 In order to better understand the significance of the present invention a brief 

discussion of the mathematical basis for the automatic selection of multiconductor 
transmission line macromodels will now be provided. 
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Referring now to the figures and to Figure 1 in particular, there is shown 
schematically a transmission line 10 as a two-port network with inputs Ii and Vi and 
outputs l 2 and V 2 . The port currents l x (s) and -I 2 (y) are related to port voltages V x (s) and 
V 2 (y) by the frequency-domain exponential solution of the Telegrapher's Equation: 









with O(s) = 


0 G + sC 








Jt + sL 0 . 



(eq. 1) 



The length of the single transmission line is d y R is the per-unit-length resistance, 
L is the per-unit-length inductance, G is the per-unit-length conductance, and C is the 
1 0 per-unit-length capacitance of the transmission line. 



15 



The matrix T(» = e ' 4(I>(5)] is the transfer matrix from the near end (ports 1) 12 to 
the far end (ports 2) 14 of the transmission line 10. For a single transmission line T(s) is 
a 2 x 2 matrix. For a multiconductor transmission line with n lines, the per-unit-length 
parameters are n x n matrices and T(s) is a 2n x 2n matrix. 



When the line is lossless, R = G = 0, and the delay per unit length is given by 
tq= 4LC and the total line delay is x = dx§. 



20 



Assume the transfer matrix of a lossy transmission line can be decomposed as 



T(s) ^€^ s)] = exp 



0 G + sC 
R + sL 0 



(eq. 2) 
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Then one can use a pure delay line to model the delay transfer matrix 
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and a purely resistive network to model the loss 



.Unfortunately, the above decomposition is 



not always valid. This is because, unlike the scalar case where e x+y = e x e y , in the matrix 
10 case e A+B ^ e A e B unless the square matrices A and B commute with each other, i.e., AB = 
BA. When two matrices do not commute the nonzero matrix [A, B] = AB - BA is called 
the commutator of A and B. 



For the single conductor line case, denote A = 



"0 
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and B = 
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C~ 
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. Then, 



15 [A,B] = 



GL-CR 0 
0 LG-RC 



Note that the above commutator expression is valid for a multiconductor line 
case. The order of matrix products GL, LG, CR, and RC is important. For single 
conductor line, the commutator is 0 if and only if GL = CR which is the same as LG = 
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RC. This is exactly Heaviside's distortionless line criteria, which proves the important 
result that the behavior of a distortionless transmission line can always be decomposed 
into that of a pure delay line and a purely resistive network. 



5 When the line is not distortionless, then it is desireable to find a criterion where 

equation 2 is approximately satisfied. To derive this criterion, first expand T(s) in a 
Taylor series to the 2nd-order in the line length d, i.e., 
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T(s) — exp 

= exp (-d[A + sB] 



0 G + sC 
R + sL 0 



(eq. 3) 



~l-d(A + sB) + d 2 (A + sB? 

2 

= l-d(A + sB) + cF(A 2 + sAB + sBA + s 2 B 2 ) 

2 

Similarly, T/ OJJ (^) and Tjeiayis) are expanded into 2nd-order Taylor series with 
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= exp (-dA) 
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(eq. 4) 



and 



Tdelay(s) = exp 

= exp (-dsB) 
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(eq. 5) 



2 
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Multiplying the expansions of T/ 0J5 (s) and T delays) and dropping terms of order 3 
or higher in d, results in 

Tioss&Tdeiayis) ~ I - d{ A + *B) + d 2 (A 2 + 2s AB + s 2 B 2 ) (eq. 6) 

2 

5 Define the error matrix 

E(s) = T/ 05 ^)T^/^) - TO). (eq.7) 

Then, based on the Taylor series expansions of equation 3 and equation 6, the 

error matrix can be approximated as 

Eto- ^rA, B1 . (eq.8) 
10 . 2 

It follows that when T(s) is approximated with the product Ti 0SS (s)T de iay(s) the 

approximation error is proportional to the matrix commutator [A, B]. This is in line with 

the distortionless case where the zero commutator guarantees that the error is also zero. 

Another interesting observation is that the error is proportional to s = jco, where co is the 

1 5 frequency of operation, and quadratic in the length of the line. In other words, line length 

has more impact on the approximation error than the frequency range of operation. 

Let now<e > 0 be a required upper bound on this approximation error. By taking 
the matrix norms of equation 8, using the 8 upper bound, and replacing s with yaw, 
20 where aw is the maximum operating frequency, the conclusion is that the inequality 

|| scP\A, B] || = ow || [A, B] || < s (eq. 9) 

2 

must be satisfied in order to guarantee that the decomposition error is below the error 
25 threshold throughout the frequency range. Note that the commutator matrix norm is 
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|| [A, B] || = max ( || GL - CR || , || LG - RC || ) 
and that for the single line case, this norm reduces to 

||[A,B]||= |LG-RC|. 

5 

Another simplification is the case when the per-unit-length conductance is zero 
(e.g., on-chip transmission lines) in which case || [A, B] || = RC and the error criterion 
becomes 

10 || ttfrA. B1 || = co^RC<s 

2 2 

Based on equation 10, define the critical line length as 

(^critical ~~ 

0) max RC) 

1 5 This is the length above which the approximation error becomes unacceptable. 

To recover the accuracy, the line must be segmented into segments of length less than the 
critical length. Some implementations of the method of characteristics require such 
segmentation especially for very lossy lines. This segmentation has a negative impact on 
the efficiency of the transmission line macromodel. An interesting observation consistent 

20 with past heuristic experience is that the critical length for the method of characteristic 
macromodel is inversely proportional to the resistive losses per unit length of the line. 
Furthermore, higher frequencies require smaller line segments. 



(eq. 10) 



(eq.ll) 
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An alternative interpretation of equation 10 is that delay extraction is needed only 
if the length of the line is above the critical length. For lines with d < d critical a rational 
approximation macromodel satisfies the required accuracy. 

Having now set forth the basis for the selection of the transmission line length at 
which either one of the delay extraction macromodel 6r the rational approximation 
macromodel is the preferred macromodel for simulating transmission lines, preferred 
embodiments of a system and method for automatically making the selection of which of 
the macromodels to use for simulating a multiconductor transmission line will now be 
described. 

Figure 2 is a schematic representation of a typical electrical cell 20 of a rational 
assumption macromodel of a transmission line. The full macromodel is a cascade of such 
cells. The cell circuit element values of Rs 22, Ls 24, Gs/2 26, and Cs/2 28 are derived 
1 5 from the transmission line length d, the number of cells and the per-unit-length R, L, G, 
and C values. The rational assumption macromodel of a transmission liife is one 
preferred approach for simulating a transmission line. 

Figure 3 is a schematic representation of a typical electrical circuit representation 
20 300 of a delay extraction macromodel of a transmission line. This macromodel relates 
the voltage-current pair at the near end (Vj, Ii) 302 of the circuit 300 to the voltage- 
current pair at the far end (V 2 , 1 2 ) 304 of the circuit 300. The values Yi 306 and Y 2 308 
of the circuit300 are admittance derived from the transmission line. The values Ji 310 
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and J 2 3 12 of the circuit 300 are controlled current sources. The extracted delay of the 
macromodel is implicit in the values Ji and J 2 . The delay extraction macromodel of a 
transmission line is another preferred approach for simulating a transmission line. 



5 Figure 4 is a graphical representation summarizing the CPU efficiency of two 

preferred categories of transmission line macromodels, i.e. delay extraction macromodel 
and rational assumption macromodel. The delay extraction macromodel is more efficient 
for simulating long, low-loss transmission lines. The rational assumption macromodel is 
more efficient for simulating short, high-loss transmission lines. The crossover point 40 
10 is not known a priori 

The problem encountered with using macromodels to simulate transmission lines 
occurs when a user is required to select the more efficient macromodel to be used for 
simulating transmission lines. The present invention provides a system and method for 
1 5 the automatic selection of the more efficient macromodel for a particular transmission 
line structure. 

The flow chart in Figure 5 shows an overall CAD process 500 for simulation of 
an electrical system having several multiconductor transmission lines. First, an input 502 
20 describing the electrical circuit to be simulated (e.g., SPICE file) is provided as an input 
to the CAD system. The CAD system may be any conventional CAD system known in 
the art. Then, for each transmission line in the inputted file, an automatic macromodel 
selection is performed 504 in accordance with the process described in conjunction with 
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Figure 6 below. After applying the automatically selected macromodel, an output 506 is 
provided which is the circuit simulation result using the cells in either Figures 2 or 3 
corresponding to the selected macromodel. 

5 A preferred embodiment of an automatic macromodel selection process 600 used 

in step 504 of Figure 5 is shown in the block diagram of Figure 6. 

Transmission line parameters, including the per-unit-length R, L, G 3 and C values, 
the transmission line length d, the maximum operating frequency aw and the model 
10 error threshold 8 are provided as inputs 602 to an automatic transmission line 
macromodel selection system 604. Using the inputted parameters, an automatic 
transmission line macromodel selection decision process 604, in accordance with the 
flow chart shown in Figure 7 and described below, is performed. 

1 5 After making the macromodel selection decision, a transmission line macromodel 

is generated 606 and provided to an electrical circuit simulator program 608. Optionally, 
macromodel information from macromodel database 610 is provided to the electrical 
circuit simulator program 608 in accordance with the selected macromodel. Electrical 
terminations of the transmission line (drivers and receivers) at the near and far ends of the 

20 transmission line and other circuitry information (e.g., reference nodes for the near and 
far ends) 612 are also provided as input to the electrical circuit simulator program 608. 
The electrical circuit simulator program 608 provides output waveforms and electrical 
circuit figures of merit related to the transmission line to be simulated 614. The output 
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waveforms and electrical circuit figures of merit are used to provide the circuit simulation 
results in step 506. 

Figure 7 is a flow chart of an automatic transmission line macromodel selection 
decision process 700, the result of which is used in block 604 of Figure 6 for selecting the 
more efficient macromodel to be used to simulate the transmission line. 

The multiconductor transmission line per-unit-length R, L, G, and C parameter 
values 702 are provided as inputs to be used to compute 704 the transmission line 
distortion per-unit-length: A = max ( || GL - CR|| , || LG - RC || ). 

The computed distortion per-unit-length A from step 704, the transmission line 

length d from block 706, and the maximum operating frequency (o max from block 708 are 

provided as inputs to compute the total distortion 710: A T = A^aw . 

2 

The computed total distortion A T from block 710 and an error threshold e from 
block 712 are compared at step 714. If A T < e, the delay extraction macromodel is 
selected 716 and used in step 604 in Figure 6. If A T > e, then the rational assumption 
macromodel is selected 718 and used in step 604 in Figure 6. 

While the delay-extraction macromodel and the rational appoximation 
macromodel are preferred simulation models, it will be apparent to those skilled in the art 
that other macromodels may be included for selection according to appropriate criteria. 
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The automatic selection of efficient transmission line macromodels in accordance 
with the teachings of the present invention may be performed on a general purpose 
computer. Preferably, the automatic selection process is used in conjunction with a CAD 
system or CAD program. The automatic selection process obviates the requirement of a 
5 user to select the most efficient transmission line macromodel for a particular 
multiconductor transmission line to be simulated. 

While there has been described and illustrated a preferred embodiment of a 
system and method for the automatic selection of transmission line macromodels, it will 
10 be apparent to those skilled in the art that modifications and variations are possible 

without deviating from the broad teachings and spirit of the present invention which shall 
be limited solely by the scope of the claims appended hereto. 
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